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Using  first-principles  theory  and  experiments,  chemical  contributions  to  surface-enhanced  Raman 
spectroscopy  for  a  well-studied  organic  molecule,  benzene  thiol,  chemisorbed  on  planar  Au(l  11)  surfaces 
are  explained  and  quantihed.  Density  functional  theory  calculations  of  the  static  Raman  tensor  demon¬ 
strate  a  strong  mode-dependent  modihcation  of  benzene  thiol  Raman  spectra  by  Au  substrates.  Raman 
active  modes  with  the  largest  enhancements  result  from  stronger  contributions  from  Au  to  their  electron- 
vibron  coupling,  as  quantihed  through  a  deformation  potential.  A  straightforward  and  general  analysis  is 
introduced  to  extract  chemical  enhancement  from  experiments  for  specihc  vibrational  modes;  measured 
values  are  in  excellent  agreement  with  our  calculations. 
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Since  its  discovery  over  three  decades  ago,  surface- 
enhanced  Raman  spectroscopy  (SERS)  has  shown  signih- 
cant  promise  for  sensing  individual  molecules  adsorbed 
near  metal  nanostructures  or  substrates  with  nanoscale 
roughness  [1-3].  In  SERS,  the  conversion  of  incident  light 
into  surface  plasmons  near  nanoscale  asperities  on  metal 
surfaces,  combined  with  chemical  and  resonant  effects,  has 
been  reported  to  yield  Raman  cross  sections  increased  by 
factors  of  up  to  10^"^,  enabling  single-molecule  detection 
[4-9].  While  the  enhancement  associated  with  surface 
plasmons  can  reach  10^  [10],  remaining  increases  and 
associated  changes  in  mode  frequencies  have  been  re¬ 
ported  to  arise  from  chemical  adsorption  [4,11-16],  as 
well  as  resonant  intramolecular  and  metal-molecule  charge 
transfer  [17-19].  Despite  a  wealth  of  prior  theoretical 
studies,  these  chemical  enhancement  (CE)  mechanisms 
are  poorly  understood  and  difficult  to  quantify  [4,12-14]. 
While  for  some  vibrational  modes  CE  is  estimated 
[4,11,12,20]  to  be  about  10-100,  there  is  currently  no  clear 
picture  for  why  certain  modes  are  enhanced  more  than 
others.  Previous  theoretical  studies  have  proposed  mecha¬ 
nisms  for  this  mode  dependence  [14-19],  but  none  of  these 
models  have  been  validated  by  experiments  or  more  rig¬ 
orous  first-principles  calculations,  leaving  the  origin  of 
dominant  chemical  contributions  to  SERS  an  open  ques¬ 
tion  [4,12-14]. 

In  this  Letter,  we  use  first-principles  calculations  and 
experiment  to  explain  chemical  contributions  to  SERS 
for  (C5H5SH)  benzene  thiol  (BT)  molecule  chemisorbed 
to  Au  surfaces.  Erom  our  density  functional  theory  (DET) 
calculations  of  static  contributions  to  the  Raman  tensor,  we 
elucidate  the  vibrational  mode  dependence  [4,12-14]  of 
the  chemical  enhancement,  explicitly  relating  modes  with 
the  largest  CE  to  those  with  the  greatest  mode-induced 
shift  of  the  molecular  frontier  orbital  energy,  as  quantified 
through  a  deformation  potential.  Relative  CE  of  BT  vibra¬ 


tional  modes  measured  at  different  probe  frequencies  agree 
quantitatively  with  our  static  calculations  for  all  binding 
sites  considered.  While  the  magnitude  of  CE  for  BT  on  Au 
is  sensitive  to  binding  motif,  its  relative  value  is  not. 

DET  calculations  are  performed  using  the  Vienna 
ab  initio  simulation  package  (VASP)  and  within  a  general¬ 
ized  gradient  approximation  [21,22].  We  model  the  BT 
adsorbate- Au  substrate  system  with  an  ordered  monolayer 
of  about  one  BT  molecule  per  nm^  bonded  to  a  flat  periodic 
Au(lll)  slab.  A  400  eV  plane- wave  cutoff  and  2X2X1 
Monkhorst-Pack  ^-point  mesh  is  used  for  calculations  in¬ 
volving  Au  slabs.  Pour  binding  geometries  are  considered, 
as  shown  in  Pig.  1(b):  fee  hollow  =  0.219  eV),  ad¬ 
atom  {E^  =  0.446  eV),  hydrogenated  adatom  {E^  = 
0.795  eV),  and  bridge  {E^  =  0.192  eV),  where  E^  is  the 
calculated  binding  energy  relative  to  a  free  Au  surface  and 
gas-phase  C5H5SH  in  the  dilute  limit.  The  gas-phase  BT 
molecule  is  simulated  in  the  same  large  supercell  as  the 
slab  geometry,  using  the  T  point. 

Static  Raman  tensors  are  constructed  mode-by-mode 
using  a  finite-differences  approach,  in  two  steps.  Pirst, 
the  dynamical  matrix  of  the  system  is  generated  by  dis¬ 
placing  each  atom  along  each  Cartesian  direction  by 
0.03  A.  Vibrational  frequencies  and  corresponding  phonon 
eigenvectors  are  obtained  by  diagonalization  of  a  truncated 
dynamical  matrix  treating  only  BT  atoms,  and  Au  atoms 
directly  bonded  with  sulfur.  Second,  we  compute  the  static 
polarizability  from  a  second-order  finite-difference  expres¬ 
sion  using  a  sawtooth  potential  with  a  gradient  of  1  mV/ A, 
and  compute  its  derivative  as  a  function  of  the  amplitude 
for  each  vibrational  eigenmode.  (See  supplemental  mate¬ 
rial  for  additional  details  [23].)  Throughout  the  Letter, 
all  modes  are  labeled  with  gas-phase  frequencies  for 
simplicity. 

In  Pig.  1(a),  we  report  Raman  cross  sections  calculated 
from  the  dominant  (nonevanescent)  component  of  the 
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FIG.  1  (color  online).  (A)  Calculated  Raman  spectra  of  BT 
(from  components  only)  for  four  different  binding  geome¬ 
tries:  adatom,  hollow,  adatom  +  H,  and  bridge.  Inset  shows  our 
calculations  of  the  orientationally-averaged  gas-phase  Raman 
spectrum  of  benzene  thiol.  All  spectra  are  broadened  with  a 
Lorentzian  of  10  cm“^  width.  (B)  Binding  geometries.  Our 
supercell  consists  of  5  atomic  layers  of  Au  stacked  along 
[111]  with  16  atoms  per  layer,  with  30  A  of  vacuum.  The  in¬ 
plane  lattice  parameters  are  kept  fixed  to  their  computed  Au  fee 
bulk  value  of  4.17  A. 

Raman  tensor,  where  z  is  the  normal  to  the  surface. 
Prominent  Raman  peaks  in  the  computed  spectra  agree 
well  with  experiment  [24].  Upon  adsorption,  BT  vibra¬ 
tional  frequencies  are  altered  on  average  by  10-20  cm“^, 
in  agreement  with  experiments  (see  [23]).  From  Fig.  1(a), 
the  presence  of  the  Au  substrate  enhances  Raman  cross 
sections  for  some  modes  more  than  for  others.  The  inten¬ 
sities  of  modes  with  larger  enhancements  exhibit  a  stronger 
dependence  on  binding  site.  For  example,  while  the  inten¬ 
sity  of  the  1100  cm“^  mode  varies  by  more  than  2  orders  of 
magnitude  across  the  different  sites  considered,  the  inten¬ 
sity  of  the  3151  cm“^  mode  remains  relatively  unaltered 
and  comparable  to  the  gas-phase  intensity. 

In  Fig.  2(a),  we  compare  the  computed  charge  density 
change  induced  by  these  two  modes,  Sp{Q^)  =  p{Qn)  ~ 
p(0),  where  is  the  vibrational  eigenvector  amplitude. 
We  find  that  the  1100  cm“^  mode  induces  significant 
charge  redistribution  within  the  Au  surface,  while  the 
3151  cm“^  mode  does  not.  This  behavior  is  entirely  con¬ 
sistent  with  mode-induced  changes  observed  in  the  DFT 
electronic  structure  [Fig.  2(b)]:  The  peak  Kohn-Sham 
HOMO  energy  £'homo  shifts  noticeably  relative  to  the 
Fermi  level  Ep  with  the  1100  cm“^  mode,  whereas  the 


FIG.  2  (color  online).  (A)  Computed  isosurfaces  (at  1.65  X 
10“^  e/}^)  of  charge  density  induced  by  freezing-in  small 
amplitudes  (0.1  A)  of  the  1100  and  3151  cm~^  modes.  For 
both  modes,  gold  atoms  are  stationary,  and  charge  rearrangement 
in  the  substrate  is  induced  by  the  molecule.  (B)  Partial  densities 
of  states  projected  on  the  molecule  near  the  Au  Fermi  energy  for 
three  different  amplitudes. 

3151  cm“^  mode  leaves  the  HOMO  peak  unchanged. 
Examining  all  modes,  we  find  that  modes  that  show  larger 
enhancements  induce  a  larger  polarization  response  in  the 
substrate  (and  larger  shift  in  E’/^-E’homo)* 

To  rationalize  our  DFT  results,  we  compare  a  model 
expression  for  the  static  polarizability  and  Raman  tensor 
with  our  more  rigorous  first-principles  calculations. 
We  consider  a  single  (dominant)  term  in  an  approximate 
single-particle  form  for  the  ground  state  electronic  polar¬ 
izability  a,  one  that  includes  just  a  lone  virtual  transition 
between  the  HOMO,  ^homo^  ^  metallic  state  at 
(/>Au*  Within  this  two-state  approximation,  an  adsorbate- 
metal  interfacial  contribution  to  the  Raman  tensor  /^bt-Au 
can  be  expressed  as  [13,14] 

_a^BT-Au^  1  _M]Mj  doj^ 

dQ,  co^  dQ,  dQ,  ^ 

where  =  (^homo  I  A'/  I  ^Au)  is  a  dipole  matrix  ele¬ 
ment  in  the  Cartesian  direction  /,  and  (o^  =  E’/^-E’homo 
is  the  difference  between  their  eigenvalues.  For  modes 
that  are  not  also  IR  active,  terms  involving  derivatives  of 
induced  dipole  moments  can  be  neglected,  and  Eq.  (1) 
reduces  to  a  term  proportional  to  —  ^  .  As  noted 

previously [13, 14],  the  co^  factor  in  the  denominator  is 
consistent  with  a  sensitivity  of  the  overall  spectral 
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enhancement  to  the  alignment  of  the  frontier  molecular 
electronic  orbital  with  Ep.  Indeed,  this  factor  can  be  used 
to  rationalize  the  binding  site  dependence  of  enhancements 
shown  in  Fig.  1(a),  and  corresponding  partial  densities 
of  states  shown  on  Fig.  3(a).  Calculated  enhancements  on 
the  adatom  site,  for  which  the  HOMO  is  close  to  Ep,  are 
significantly  stronger  than  on  other  sites,  where  this  level  is 
broadened  and  further  away  from  Ep.  However,  this 
factor  is  the  same  for  all  modes,  and  cannot  explain  the 
mode-dependence  of  CE  seen  in  Fig.  1(a). 

The  strong  modification  of  the  BT  Raman  spectra  by 
the  Au  substrate  can  be  explained  via  computation  of  a 
deformation  potential,  for  each  mode,  i.e.,  the 

change  in  molecular  electronic  level  alignment,  relative 
to  Ep,  induced  by  a  particular  vibration  mode.  In  Fig.  3(b), 
we  plot  the  difference  —  /^Adatom+H  against 

Because  the  Kohn-Sham  HOMO  level  of  the 
adatom  +  H  geometry  is  much  further  from  Ep  than  for 
the  adatom  binding  site  [see  Fig.  3(a)],  is  large,  and  the 
adatom  +  H  geometry  has  a  negligible  interfacial  contri¬ 
bution.  Thus,  to  an  excellent  approximation,  the  difference 
y^Adatom  _  ^Adatom+H  j-gmoves  intramolecular  contributions 
to  the  Raman  tensor  unrelated  to  Eq.  (1).  Indeed,  in 
Eig.  3(b),  we  find  a  remarkable  correlation  between  the 
interface  contribution  to  the  deformation  potential  and 
enhancement  for  most  of  the  vibrational  modes.  The 
modes  that  deviate  from  the  linear  trend  at  407  cm“^ 
(Au-S  stretch),  1475  cm“^  (phenyl  ring  stretch),  and 
3151  cm“^  (C-H  stretch)  have  significant  IR  activity. 


FIG.  3  (color  online),  (a)  Partial  densities  of  states  projected 
on  the  molecule  for  different  binding  sites  shown  about  the  Au 
Fermi  energy,  (b)  Difference  /^Adatom  _  ^Adatom+n  versus 

dconldQn.  (c)  The  deviation  8  of  these  modes  from  linear  trend 
in  (b)  plotted  versus  their  contribution  to  the  lattice  susceptibil¬ 
ity,  which  is  proportional  to  the  square  of  the  ratio  of  the  mode 
dynamical  charge  and  its  frequency. 


In  the  static  limit  used  here,  polarization  induced  by  these 
IR-active  modes  screens  the  electric  field  experienced  by 
the  molecule,  leading  to  a  reduction  in  their  Raman  cross 
sections  [25]  by  a  factor  proportional  to  (Z*)^/ col,  with  the 
mode  dynamical  charge  given  by  Z*  =  Ayit/AQ^  and 
where  Ayu.  is  the  mode-induced  interfacial  dipole  moment. 
In  Fig.  3(c),  we  show  the  deviation  S  of  the  IR-active 
“outlier”  modes  from  the  linear  trend  observed  for  non- 
IR-active  modes  versus  (Z*)^/ co^.  The  correlation  between 
S  and  (Z*)^/ co^  confirms  that  modes  with  larger  contribu¬ 
tions  to  the  screening  lead  to  greater  deviations  from  the 
two-state  model  [Fig.  3(b)].  We  note  however,  that 
although  important  in  our  static  calculations,  this  screening 
effect  will  have  impact  only  below  infrared  frequencies, 
and  thus  will  be  inconsequential  for  typical  probe  frequen¬ 
cies,  where  the  local  fields  will  vary  too  rapidly  for  the 
IR-active  vibrations  to  respond. 

With  this  information,  we  can  now  connect  the  strong 
modification  of  Raman  spectra  by  substrates  to  mode- 
specific  changes  of  the  electronic  structure  of  the  metal- 
adsorbate  interface:  modes  with  the  largest  interfacial 
contribution  to  the  change  in  polarizability,  as  quantified 
through  a  deformation  potential,  result  in  the  most  sub¬ 
stantial  chemical  enhancements.  For  BT,  these  modes  are 
those  that  break  the  conjugation  of  the  HOMO  (supple¬ 
mental  material  [23]).  Figure  2(a)  shows  vividly  how 
the  1100  cm“^  mode  breaks  the  resonant  character  of 
the  carbon  ring,  while  the  3151  cm“^  mode  leaves  the  it 
symmetry  of  the  electrons  on  the  phenyl  ring  intact.  This 
suggests  that  for  future  adsorbates,  the  nature  of  modes 
with  the  largest  CE  might  be  intuitively  rationalized  a 
priori. 

To  validate  our  theory  of  chemical  enhancement, 
we  compare  with  experimental  Raman  and  SERS 


Calculated  relative  enhancement 


FIG.  4  (color  online),  (a)  Au  SERS  substrate  used  in  our 
measurements  with  BT.  (b)  Raman  and  SERS  BT  spectra, 
(c)  Calculated  relative  enhancements  versus  experimental  rela¬ 
tive  enhancements  for  adatom,  bridge,  and  hollow  sites. 
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measurements  of  BT  on  rough  Au  substrates  [Figs.  4(a)  and 
4(b)].  SERS  substrates,  consisting  of  roughened  SiGe  sur¬ 
faces  [Fig.  4(a)]  coated  with  30  nm  of  Au,  were  incubated 
in  3  mM  BT  solution  (Sigma  Aldrich  W361607)  in  metha¬ 
nol  overnight,  then  gently  rinsed  with  methanol  and  dried 
by  nitrogen  gas.  Raman  (neat  solution)  and  SERS  spectra 
were  collected  at  two  wavelengths,  632.8  nm  and  785  nm, 
using  an  inverted  microscope  setup  coupled  to  a  spec¬ 
trometer  (Acton  SpectraPro  2300i)  equipped  with  a 
liquid-nitrogen-cooled  charge-coupled  device  camera. 
To  eliminate  uncertainties  associated  with  the  number  of 
molecular  analytes  in  comparing  Raman  and  SERS  inten¬ 
sities,  we  normalize  ratios  of  Raman  and  SERS  spectra 
peak  heights  to  the  996  cm“^  mode,  which  has  only  a 
modest  enhancement  from  our  calculations.  We  note  that 
absolute  CEs  can  be  obtained  if  we  normalize  to  a  mode 
with  zero  deformation  potential,  dconldQ^,  e.g.,  mode 
1348  cm“^  in  Fig.  3(b).  As  the  1348  cm“^  mode  is  not 
easily  observed  experimentally,  we  use  996  cm“^  in  this 
case.  Assuming  the  electromagnetic  enhancement  is  the 
same  for  all  modes,  this  relative  enhancement  will  reflect 
CE.  This  assumption  is  acceptable  for  modes  within  a  few 
hundred  wave  numbers  of  the  996  cm“^  mode,  based  on 
the  width  on  relatively  low  Q  factor  for  localized  plasmon 
resonance  in  Au  or  Ag.  In  a  more  general  case,  the  elec¬ 
tromagnetic  enhancement  is  not  the  same  for  all  vibrational 
modes,  and  a  plasmon  resonance  dispersion  correction 
needs  to  be  applied  [26].  Additionally,  in  order  to  compare 
data  from  different  experimental  setups  or  different  exci¬ 
tation  wavelengths,  spectral  responses  of  specific  experi¬ 
mental  setups  must  be  taken  accounted  for  (see  [23]). 

In  Fig.  4(c),  we  compare  Raman  cross  sections  from 
Fig.  1(a)  (computed  from  to  averaged  solution  phase 
data  from  Fig.  1(a)  (inset).  Binding  geometries  for  BT  on 
flat  Au  surfaces  are  the  subject  of  debate  in  the  literature; 
and  for  rough  surfaces  at  room  temperature,  a  variety  of 
adsorption  sites  will  be  available,  and  subject  to  thermal 
fluctuations  [9,27].  By  comparing  three  very  different  op¬ 
timized  binding  sites  to  the  experiment,  we  sample  differ¬ 
ent  possibilities  for  time-averaged  experimental  binding 
morphologies.  Surprisingly,  all  three  geometries  show 
good  correlation,  with  the  adatom  site  perhaps  showing 
the  most  linear  trend.  (We  note  that  SERS  measurements 
do  not  observe  the  S-H  stretching  mode  at  2600  cm“^ 
indicating  the  loss  of  hydrogen  at  the  S-Au  in  the  room- 
temperature  experiments.)  Importantly,  the  relative  theo¬ 
retical  enhancements,  also  normalized  to  the  996  cm“^ 
mode,  are  in  excellent  quantitative  agreement  with  the 
present  experiments,  as  well  as  others  in  literature  taken 
from  different  substrates  [28,29]  (supplementary  material), 
indicating  standard  DFT  approaches  can  capture  accurately 
and  quantitatively  dominant  contributions  to  CE,  even 
in  the  static  limit.  The  quantitative  agreement  between 
theory  and  experiment,  and  the  connection  between  defor¬ 
mation  potential  and  chemical  enhancement,  provides  new 


opportunities  for  detection  and  control  of  adsorbate-metal 
interactions  through  SERS. 
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